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<D . Abstract 

We have found the analytical solution of the LO-evolution equation for off-forward distributions 
^ ■ which arise in the processes of deeply virtual Compton scattering or exclusive production of 
mesons. We present the predictions for their evolution with an input distribution taken from 
recent bag model calculations. 



1. Introduction. Recently, X. Ji proposed to explore deeply virtual Compton scattering (DVCS) 
(see Fig. |]) to get deeper insight into the spin structure of the polarized nucleon and, namely, 
to determine the total quark angular momentum (spin and orbital angular momentum) in the 
nucleon [|l|, [|. It turned out that this process is very interesting in its own right, since it allows 
to obtain information about some non-perturbative off- forward parton distributions (OFPDs), 
which are inaccessible in ordinary inclusive measurements such as deep inelastic scattering (DIS). 
Planning of experiments to measure DVCS is under way, but it is hampered by the lack of reliable 
theoretical predictions in the interesting (u, r], Q 2 )-range. A crucial step towards this goal is to 
understand the Q 2 -evolution of OFPD. In the present short note we, therefore, present the Q 2 - 
evolution of the nonpolarized quark distribution function in the leading order (LO) approximation 
starting form MIT bag model predictions at a low energy scale ||. 
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Figure 1: Generic form of the DVCS amplitude. 

To make this paper self-contained, as much as possible, let us first recall the main points of 
the perturbative (LO) analysis ||, [1], |], |7|, || ^, [1(J of DVCS amplitude. In the configuration 
space such a Compton process, given by the following expression 

T^(u, V ,Q 2 )=zJ d*ze* zQ (p 2 \T {J» (-z/2) J» (z/2)}\ Pl ), (1) 

is dominated by singularities of the T-product of electromagnetic currents J M on the light-cone. 
In momentum space this dominance corresponds to the generalized Bjorken limit: 

— Q 2 — ► oo, PQ — > oo, u, A 2 -fixed, 

where the variables are defined as follows Q — ^(qi + Q2), P — (Pi + P2), A = p\ — p 2 = q 2 — Qi- 
DVCS is characterized by two scaling variables u = — ^3 and r] = Q, the former is analogous 
to the Bjorken one in ordinary DIS while the latter gives the magnitude of the skewedness of the 
process in question and parametrizes different high-energy two-photon reactions according to the 
Table 0. 

Notice that in general the unpolarized tensor T^ u can be decomposed into four (three in the 
case of DVCS) independent kinematical structures that separately respect current conservation 
(similarly for polarized hadrons) |TT| . However, only two of them survive the general Bjorken limit, 
while the remaining ones are kinematically suppressed and contain higher twist contributions. In 
LO only one kinematical structure is independent, so that the longitudinal helicity amplitudes 
vanish. This is analogous to the Callan-Gross relation in unpolarized DIS. Motivated by the 
result that only transverse helicity amplitudes contribute in LO, we can write, e.g. for unpolarized 
scattering 

T(A',A) = e^(A / )T^(A)^^(A')er(A)5 A1/2 5v 1/2 T w . (2) 
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Table 1: The value of the scaling variable rj for different two-photon processes 
Adding the surviving spin-dependent structure, we getQ: 



ut-l + iO ujt + 1 - zO J ' 

IrS = %r T^E^a^A 2 )! }. (3) 

[H 2PQ J-i Y V 'W-l+iO cui + l-zOj v; 

Here the off-forward distributions are defined in terms of the light-cone Fourier transformation of 
the renormalized twist- 2 light-ray operators: 

| ^ | (t,r],A 2 ,fi 2 ^ = J ^ e i" p +<p 2 M-K/2n) j ^ | V(«/2n) |^ |pi)a + =,p + , (4) 

and we have omitted the path ordered link factor which ensures gauge invariance. The integration 
range of the variable t in Eq. (Q) results from the spectral properties of the off-forward distribution 
amplitude that can be derived in a straightforward manner with the help of Jaffe's approach JT3 



for studying the support properties of parton densities. In our conventions the variable t has 
no direct interpretation as longitudinal momentum fraction, but it has the advantage that its 
region is process independent: — 1 < t < 1. The same feature is appropriate to the quark 
^(x) (0 < x < 1) and anti-quark JF"(x) (—1 < x < 0) non-forward distributions introduced by 
Radyushkin []5j. The (^-dependent limits appear only when one attempts to combine them into 
one function F?(x) = 9{x)6{x)!Fq{x) — 9(( — x)9(x + 1 — C)^?(C — x ) which is exactly the OFPD 
introduced above in Eq. in different notations. The non-diagonal distributionF] f g (xi,X2) 
of Collins, Frankfurt and Strikman is parametrized in terms of the momentum fraction of the 
considered partons and consequently their limits depend on the skewedness of the processes. Thus, 
the momentum fractions of the incoming and outgoing partons are x = and x — £ = 

respectively, with £ = For t > r\ both momentum fractions are positive and, therefore, 

a quark field with longitudinal momentum xpi arise from the incoming hadron and enters with 
momentum (x — Qpi into the outgoing hadron. For — T) > t both momentum fractions are negative 
and the analogous picture for two anti-quarks holds true. In these regions, the corresponding off- 
forward quark and anti-quark distributions merge in the forward case to the well-known quark and 
anti-quark distributions with support < t < 1 and — 1 < t < 0, respectively. For —rj < t < 77, 
x is positive while x — ( is negative, so that both partons come from the same hadron. Here the 
OFPD looks like a distribution amplitude. Indeed, inserting the vacuum and one-meson state 

1 Recently, all NLO coefficient functions |l2|, [l(J as well as a s -corrections to the eigenfunctions of the two-loop 
QQ- and QG-evolution kernels became available. 

2 In the Radyushkin's conventions f 9 (x, x — £) = (x)/x(x — C) H ■ 
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in the definition (|4]) instead of incoming and outgoing nucleon state, we recover the definition of 
the meson distribution amplitude depending on the momentum fraction < x — (1 + i)/2 < 1, 



where r] = 1 [Q. Thus, OFPDs are hybrids and their probability interpretation depends upon 
the interplay between the respective values for t and r\. 

Finally, to end up with this introduction, we point out a few difficulties of the measurement 
of the OFPDs in DVCS. DVCS can be explored in the reaction 



(5) 



However, DVCS and Bethe-Heitler (BH) sub-processes have the same final states. Hence, they 
interfere and the measured cross section contains the squares of pure DVCS and BH parts as well 
as the interference term: 



da oc |T DV cs| 2 + I^bh| 2 + (^dvcs^bh + ^bh^dvcs) 



(6) 



Fortunately, the contribution from the BH process is known from the measurements of the nucleon 
form factors. In unpolarized or double spin experiments DVCS amplitude is given as convolution 
[compare with Eqs. (|2|) and (|3])], so that a direct extraction of the OFPDs seems to be very 
difficult^. In such experiments the BH contribution dominates over DVCS for small A 2 . Its 
importance decrease, however, with increasing A 2 [|lj. In single spin experiments a non- vanishing 
cross section arises from the imaginary part of the amplitude. Since BH has only a real part, the 
single spin asymmetry is proportional to the imaginary part of the DVCS amplitude, so that the 
OFPDs can be measured directly at the point t = ±l/a>. Hopefully, as a first step in the direction 
of DVCS measurements, the predicted scaling behaviour in the generalized Bjorken regime will be 
investigated at CEBAF gj, |. 

2. Evolution equation. The evolution of OFPDs arises technically from the renormalization 



procedure for the light-ray operators [15], which satisfy in the non-singlet channel the well-known 
renormalization group equation [O, 17, 18, [HJ. Employing the definition (01), it is straightforward 



to derive from the latter the evolution equation for the OFPD: 



dQ 



1 dt' 




g NS (t' )?7; A 2 ,Q 2 ). 



(?) 



In the restricted region \t,t'\ < 1 the kernel V(t,t') coincides with the Efremov-Radyushkin 

V{2x-l,2y-l 



Brodsky-Lepage (ER-BL) one 



'0<x,y<l • 



(8) 



Moreover, the continuation to the whole t, t'-plane is unique [§], so that the extended evolution 
kernel V(t, t') can be easily recovered from the corresponding ER-BL analogues. From the support 
restrictions and charge conjugation symmetry a representation was derived in p0| , which is valid 
in leading order: 

V(t,t') = Q° u (t-t',t-l)V(t,t') + (t^-t,t' (9) 

where 05 1 (xi,x 2 ) = (9(xi)9(—x 2 ) — 9(x 2 )0(—Xi)) / (x\ —x 2 ), and V{t,t') is an analytic function of 
t and t'. Thus, the extended kernel can be obtained from the exclusive ER-BL one by replacing 

3 The imaginary part of the LO scattering amplitude (||) is expressed in terms of the OFPDs at t — ±l/u>, while 
the real part is given as a Hilbert transformation and can be inverted in principle. 
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the ordinary ^-functions in front of V(t,t') by the generalized ones, namely, 8(t — t')/(l — t') 



©11 (t — t',t— 1). For concrete examples the reader is referred to Ref. |T0|. Beyond LO, further 
contributions appear in Eq. (^|), but they are not of relevance for our present consideration. 

For r] — > the evolution equation ([?[) turns into the DGLAP one with a splitting function 
given by the limit: 



P(z) = lim 
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(10) 



3. Solution of the evolution equation. The evolution equation analogous to Eq. (|7]), but for gluons, 
was attempted to be integrated numerically ||. However, it remains an important issue to find an 
analytical solution. The way to do that is based on the use of the conformal invariance of the free 
field theory which ensures that conformal 2-particle operators do not mix under renormalization 
at LO jTJ], |2T], [22|, |TTJ. Such conformal operators are given in terms of Gegenbauer polynomials 



moments 




ip and their off-forward matrix elements can be obtained from the 



A 2 , Q 2 ) 



rf 



dtCf (-N(i,77,A 2 ,Q 



P: 



j+ 



i W*>ii(0) \ Pl ). 



In the following, we restrict ourselves to the flavour non-singlet channel since up to now no 
estimation of the off-forward gluon distribution has been performed yet as long as available strong 
interaction models do not contain the gluon fields at all. The generalization of our results to the 
singlet channel is straightforward. The only complication which is due to quark-gluon mixing can 
be treated in the same way as in DIS. For r] = 1 it is well-known that the Gegenbauer polynomials 
diagonalize the ER-BL-kernel. Moreover, the support of the extended evolution kernel ensures 
this property for arbitrary 77: 






(12) 



The remaining problem is to find the inverse Mellin transformation of the moments (0). The 
support property^ of g NS (t, 77, Q 2 ) allows to expand it with respect to an appropriate complete set 
of polynomials C%(t) which are orthogonal in the domain —1 < t < 1 with the weight function 
w{t\v) = (1-t 2 ^- 1 - 



1 

' 2 ■ 



q" S (t, V ,Q 2 )=f:^c}(t)a J (r ] ,Q 

j=0 1S j\2> 



213) 
2 I ! 



(13) 



2 2u+1 wrir is a normalization factor. It is straightforward to calculate the 
expansion coefficients in terms of the conformal moments ([TT|): 



where Nj{y 



a j(v,Q 2 \l) = Y, a Av\li)<lk S (v,Q 2 o)ew\lk f 



2 2 da ct s (cr) 
g V 2tt ^ 



(14) 



4 In the subsequent discussion we omit, for brevity, the dependence of the distributions on the momentum 
transferred squared. 
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where jj = C F [—2ip(j + 2) + 2-0(1) + l/(j + 1) - l/(j + 2) + 3/2] are the well-known forward 
non-singlet anomalous dimensions, and %fe(^|f) can be written in a very compact manner by 
employing the definition of hypergeometric functions^]: 



a jk (rj\u) 



1 Azm^imw 



1 T] k N k {u) 



7 Jk 



i + (-iy- k 



(-i)^r(i/ + i±* 

i> + fc)r(i + ^) Zii V ^ + ^ + i 



2 i ~ 2 



(15) 



Let us turn to the consideration of the limiting cases. Since 2F\{— n, a + n; a + 1; 1) = <5„o, we 
obtain for rj — 1: 



«fc(i,Q 2 ll) 



exp < 7 fc 



2 2 dcr a s (cr) 

'2 cr 27T 



(16) 



and, therefore, from Eq. ( |T3"D we recover the well-known LO result for the evolution of the meson 
distribution amplitude [Q. In the limit rj — > 0, the solution of the DGLAP equation expanded 

3 

with respect to Gegenbauer polynomials C fc 2 , has the following coefficients 



%(0,Q 2 



E 



4 )^r(2±p) 



-4-gf (0,gyexp<| 7fe 



^ r(3^)r(2±|^) 

A; — j even \ / \ / 



2 da a s (a) 
Ql a Air 



(17) 



In this limit conformal moments coincide up to an overall normalization with the usual ones in 
DIS: 



2 k T(u + k) f l 
k\T{u) 



g NS (x,0,Qg) + (-l)V S {~x,0,Ql 



For asymptotically large Q 2 all conformal moments with k > will be suppressed due to non- 
zero anomalous dimensions and, consequently, only j-even expansion coefficients survive in this 
limit: 



lim aj(r],Q\l 

(-i)^r(^) 
r(l)r(¥ 



(19) 



-j/2,(3 + j)/2 
5/2 



?7 Uo l^Qo) for 3 even- 



It is not hard to check that they belong to the following asymptotic off-forward distribution which 
was originally found in Ref. ||: 



(20) 



It turns into <5(t) when 77 — > 0. From the asymptotic form found so far, we can argue that, in 
general, OFPDs will be enhanced in the region < \t\ < rj and suppressed for 77 < \t\ < 1 when 
evolved upwards in Q 2 . To demonstrate these effects in a clear manner, we choose a ^-independent 
toy-model input distribution 3/4(1— t 2 ) (which is in fact unphysical for arbitrary 77). For 77 = 1 this 
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Figure 2: The toy input OFPD (solid line) given at Q 2 = 0.5 GeV 2 is evolved up to Q 2 = 5 GeV 2 , 
where A QCD = 220 MeV. The dashed and dotted lines show the evolution for rj = 0.5 and rj = 0, 
respectively, while for rj = 1 the input distribution does not evolve since it has already the 
asymptotic shape. 



distribution possesses already the asymptotic form, while for smaller r/'s the evolution provides the 
mentioned behaviour as illustrated in Fig. ^. The evolution has been done with the help of Eqs. 
(|I1I|)-(|15|), while the conformal moments of our toy distribution at the input scale Q 2 = 0.5 GeV 2 
was computed from Eq. (|TTJ) . Taking the first 80 terms in Eq. flT3|) the OFPD was evolved to 
Q 2 = 5 GeV 2 for r] = {5/10,0} and t = 0, ±1/20, ±2/20, 1. To avoid numerical problems 
arising from the oscillation of the Gegenbauer polynomials the calculation were done exactly; and, 
finally, the result was interpolated to a smooth function. 

4- Q 2 -evolution of the bag model motivated OFPD. A first non-perturbative estimation of the 
off-forward valence quark distribution function has been done using the MIT bag modelQ ||. In 
this model calculation it turned out that at a scale /x^ag — 0-2 GeV 2 the ^-dependence of the 
off-forward u and d quark distributions is extremely weakQ and these functions vanish for negative 
t as well as for t — > 1. (The nonzero result of the calculated t-dependence in the vicinity of unity 
is a model artifact, which reflects the fact that the incoming and outgoing protons are not good 
momentum eigenstates.) For positive t the distributions are positive with a maximum at t ~ 0.4. 
To be able to satisfy the sum rules which are obeyed by the off-forward distributions, namely, the 
first two moments should be independent of the skewedness parameter 77 

Q(tiV,Q 2 ) = Jjfol { q^[o) }bi>. (21) 

(Qqw j s a quark p ar t of the energy-momentum tensor), we neglect the ^-dependence altogether, 
which is not far from the bag model results. Since the main goal of the present study is to 
acquire some intuition about the evolution properties of OFPD, we do not pursue the aim for 
construction of more realistic models; rather motivated by the given results, we take a very simple 
parametrization of the OFPD in the non-singlet channel, namely 

q NS (t, ?7, Ql = ) = 60 9{t)t 2 {\ - t) 3 , (22) 




5 To be general we have not specified the index v. for the quark operators v — 3/2, for gluons v = 5/2. 
6 Quite recently there appeared a paper 24 with evaluation of the OFPDs in the instanton inspired chiral 



quark-soliton model of the nucleon. However, it is not possible to use their results here as only flavour singlet 
quark density was considered. Moreover, flavour non-singlet combination is l/Af c -suppressed in this approach. 
7 This is in contrast with the result of Ref. p4j, where the ^-dependence is rather strong. 
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with a first moment normalized to unity. 

In Fig. we evolve this input as described previously (we took also discrete values for rj = 
0, 1/10, . . . , 1 and interpolated the result with respect to t and rf) up to the scales Q 2 = 2 GeV 2 , 
Q 2 = 200 GeV 2 , and asymptotically large Q 2 ; for other parameters, we set Nf — 3, Aqcd = 
220 MeV in a s (Q 2 ). As expected from our previous discussion, the distribution spreads in t 
for larger value of rj and shrinks for smaller ones. These graphs clearly show that during the 




Figure 3: The evolution of the bag model motivated OFPD (|22"D for Aqcd = 220 MeV, where the 
scales are Ql = fil ag (a), Q 2 = 2 GeV 2 (b), Q 2 = 200 GeV 2 (c), and asymptotically large Q 2 (d). 

evolution the large-?] part spreads over the whole range of the momentum fraction while the 
small-// distribution is pushed towards t — *■ and concentrates in the vicinity of zero. Since we 
have omitted in the present analysis the gluon sector, we can expect that taking it into account we 
will get, probably, an even more enhanced function in the small-t (and rj — > 0) region for moderate 
values of Q 2 . As is easily seen from the above figures with growing scale the OFPD approaches 
the asymptotic form. Namely, for r\ — > 1 it takes the (1 — t 2 )-shape of the asymptotic distribution 
amplitude, while for small 77 it becomes 5(t). For intermediate rj it smoothly interpolates between 
these limits. This is a quite general feature that should be obeyed by any reasonable (read physical) 
off-forward parton density. 
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